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Abstract 

The specific heat of electron-doped Pri.gsCeo.isCuC^ single crystals is re- 
ported for the temperature range 2 - 10 K and magnetic field range - 10 T. 
A non-linear magnetic field dependence is observed for the field range - 2 T. 
Our data supports a model with lines of nodes in the gap function of these 
superconductors. Theoretical calculations of the electronic specific heat for 
dirty d-wave, clean d-wave, and s-wave symmetries are compared to our data. 
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The order parameter (gap) symmetry of the high-Tc cuprate superconductors(HTSC) 
is an important parameter in attempting to understand the pairing mechanism in these 
materials. For hole-doped cuprates experimental evidence strongly favors <i-wave symmetry 
]1].|2|]. Surprisingly, early experiments on electron-doped(n-type) Ndi.ssCeo.isCuC^ (NCCO) 
suggested a s-wave symmetry. Recent penetration depth ]3[], tri-crystal ||, photoemis- 
sion j5|, Raman scattering || and point contact tunneling experiments @ on NCCO and 
Pr2- a; Ce a ;Cu04 (PCCO) favor a <i-wave symmetry. In addition to these measurements which 
show s-wave or d-wave symmetry, there are penetration depth || and point contact tun- 
nelling experiments that have shown evidence for a change in the order parameter as 
the doping changes from under-doped(d-wave) to over-doped (s-wave). However, since these 
prior measurements on the n-type cuprates are surface sensitive there is a need for bulk 
measurements (e.g. specific heat) to convincingly determine the pairing symmetry, as was 
the case for the p-type cuprates P-|T2|]. 

The specific heat is sensitive to low temperature electronic excitations. Different gap 
symmetries have different density of electronic states close to the Fermi level. Conven- 
tional low-T c superconductors show a s-wave gap symmetry in which the electronic specific 
heat has an exponential temperature dependence, C e i oc T 1 - 5 e~ A ' kT , where A is the energy 



gap For a clean (i-wave superconductor electronic excitations exist even at the lowest 
temperatures. The electronic DOS is predicted to have a linear energy dependence close to 
Fermi level, and this shows up in the electronic specific heat as C e \ oc T 2 



In the mixed state, there are two types of quasiparticle excitations in the bulk of the 
superconductor: bound states inside the vortex cores, and extended states outside the vor- 
tex cores. In conventional s-wave superconductors, the in-core bound states dominate the 
quasiparticle excitations, therefore the electronic specific heat is proportional to the number 
of vortices. The number of vortices is linear in field, therefore the electronic specific heat is 
also linear in field. ||15|| . In a superconductor with lines of nodes(e.g. ci-wave symmetry), 
the extended quasiparticles dominate the excitation spectrum in the clean limit. It has 



been shown that the electronic specific heat has a yH dependence in the clean limit [16 



at T=0. For non-zero temperatures there is a minimum field that depends on temperature 
after which the y/~H dependence should be observed. In the dirty limit the energy scale 
related to impurity bandwidth (or impurity scattering rate) is much larger than the energy 
scale related to the Doppler shift due to magnetic field(the dominant mechanism for the 
clean d-wave case), and much less than the superconducting gap maximum. In this limit, 
i.e. ksT « (H/H c2 )A « hjo « A , where A is the gap maximum and Hjq is the 
impurity band width, the magnetic field dependence deviates from V~H, and an H\og(H) 
like dependence is predicted below a certain field H*, which depends on temperature and 



impurity concentration in the sample [17 



In this paper we present the first magnetic field dependent specific heat measurements 
on n-type cuprates which probe the symmetry of the superconducting state. The electronic 
specific heat has been observed to have a non-linear magnetic field dependence. The theo- 
retical model for a clean <i-wave symmetry fits reasonably well to our data, however there 
are deviations from this type of field dependence below H*=0.G T (Fig.3). We find that a 
H\og(H) type dependence gives a better fit to our data over the whole range, which means 
our data can better be described by a dirty d-wave symmetry. It is important to emphasize 
that the main point of this work is to address the question of s-wave vs d-wave, rather than 
clean d-wave vs dirty d-wave. 

The specific heat data was obtained in the temperature range 2 - 10 K and the magnetic 
field range - 10 T using the relaxation method [|ll|. The measurements were repeated in 



two systems, a home-made setup and a Quantum Design PPMS with some modifications 
on the sample holder to remove the field dependence of the original chip. The addenda 
consists of a sapphire substrate with a thermometer and heater, and Wakefield thermal 
compound to hold the PCCO crystal. The addenda was measured separately and was found 
to have no magnetic field dependence within the resolution of our experiment(±2.5%). The 
experiment was done on several optimally doped Pri.ssCeo.isCuC^ single crystals (the mass 
of the crystals was 3-5 mg). The sample heat capacity is approximately equal to two times 
that of the addenda at T = 2 K, and equal to that of the addenda at T = 10 K. The crystals 



were grown by the directional solidification technique [[HJ. The samples were characterized 
with a SQUID magnetometer and found to be fully superconducting, with similar transition 
temperatures T c =22 K±2 K. 

The specific heat of a d-wave superconductor usually has the following main contribu- 
tions: the electronic contribution, which could have the form 7T or 7T 2 depending on the 
field and temperature range the measurement is done, the phonon contribution, which at the 
temperature range of our experiment can be written as /3T 3 , and a Schottky contribution, 
which is caused by spin-1/2 paramagnetic impurities ||20|| . Furthermore 7 = 7(0) +'~f(H), 
where ^{H) gives the field dependent part of the electronic specific heat coefficient, and 
7(0)T is the residual linear temperature dependent part of the electronic specific heat. 7(0) 
is sample dependent, and its origin is not completely understood. f^-pH- Non-electronic 
two-level systems away from the copper-oxide planes are one of the possible candidates for the 



origin of this term [ 17| . This term has been found in all hole-doped samples studied PHTT 



Fig.l shows temperature dependence of the specific of PCCO heat at four different fields, 
T, 1 T, 2 T, and 10 T applied perpendicular to the ab-plane of the crystal. The field range 
0-2 tesla is the relevant field range to extract the gap symmetry information ]TB |, and at 



H=10 T the sample is completely in the normal state(if C 2 = 8T at T=2K). Driving the 
sample to the normal state enables us to extract an important parameter, 7„ = 6.7 ± 0.5 
mJ/mole-K 2 , which is needed to compare our data to theoretical predictions quantitatively. 
A global fit which assumes the phonon coefficient, j3, constant for all fields and 7 variable 
gives a 7(0) = 1.4 ± 0.2 mJ/mole-K 2 . This value of 7(0) is consistent with the values found 
for 7(0) in the hole-doped superconductors(7(0) ~ 1 mJ/mole-K 2 for YBCO |p Hil|l ). The 
fact that we do not have any Schottky upturn at low temperatures for any field shows that 
our sample is free from a detectable level of magnetic impurities. From the slope of the 
lines, obtained through a global fit, (3 = 0.29 ± 0.01 mJ/mole K 4 , and a Debye temperature 
Od = 362 ±4 K has been extracted. These values are in reasonable agreement with the 
other published data in the literature(/3 = 0.244 mJ/mole K 4 , and 6d = 382 K J2~T||). 

Since the phonon specific heat is field independent and there is no Schottky contribution 



to the specific heat, subtracting the zero field specific heat from the specific heat at other 
fields gives the field dependent part of the electronic specific heat. Fig. 2 shows the field 
dependent part of the electronic specific heat,7(if )T, vs magnetic field at 3.4 K in the field 
range - 8 T. Fig. 3 shows theoretical fits to the 3.4 K data in the field range - 2 T. The 
clean d-wave fit is calculated using the equation: ITT] 



c -=^(i)(^) 1/, ^(^)< 1 « 

where 7„ = 6.7 mJ/mole K 2 (from the intercept of our 10 T data in Fig. 1), H C 2 = 10 T, 
and a = 0.7(this value was found experimentally for YBCO ||) are used. The clean d-wave 
fit is clearly better than the linear s-wave fit. Possible non-linear behavior in an s-wave 
superconductor will be discussed below. 

Even though the clean d-wave fit has a much better consistency than the linear s-wave fit, 
there seem to be deviations between our data and the clean d-wave theory below H* = 0.6 T, 
which would be expected from a dirty d-wave superconductor. In fact, being in the dirty 
limit is not unexpected, since the penetration depth measurements performed on similar 
crystals, grown by this group, were also consistent with dirty d-wave symmetry ||. For 
clean d-wave symmetry the change in the penetration depth as a function of temperature 
is linear in temperature(AA(T) oc T), whereas a quadratic temperature dependence(AA(T) 
oc T 2 ) is expected if the nodes are filled by impurity states, i.e. dirty d-wave. A quadratic 
temperature dependence has been observed consistently by two different groups on many 
crystals they studied 0. If a dirty d-wave function of the type C e i(H) = Alog(B/H) is 
fitted to our data, an excellent fit is obtained for the fitting parameters A = 6.2 ± 0.6 
mJ/mole-K and B = 17.6 ± 4.7 tesla. We compared these parameters with the theoretical 
predictions ]17j calculated from the equation: 



where Tr/o is the impurity band width, A is the superconducting gap maximum, and a is a 
geometrical factor related to vortex lattice geometry. Substituting H C 2 = 10 T and a = 1, 



the fitting parameter B = ( 2 ^ 2 -) can be calculated to be 15.7 tesla, in reasonably good 
agreement with the value generated by our data B = 17.6 ±4.7 tesla. The other fitting 
coefficient A can be calculated from A = InT ( jj I /a 2 )- However the parameter 70, 
which is related to the density of impurities in the sample and is sample dependent, is not 
known. Therefore by using the experimental A = 6.2 ± 0.7 mJ/mole-K value, hjo can be 
estimated to be 2.1 K, which is in good agreement with high quality YBCO samples P3J . 

Our analysis at the other temperatures(Fig.4) also produced results similar to the 
T=3.4 K data. For T=2.3 K, A = 3.1 ± 0.2 mJ/mole-K and B = 18.0 ± 4.7 tesla, and 
for T=2.7 K A = 4.3 ± 0.5 mJ/mole-K and B = 18.7 ± 4.9 tesla are found. Theoretically 
the coefficient B should be the same for all temperatures, and the coefficient A should be 
linearly proportional to the temperature. The best fits to the data generated the same values 
for B within the error range, and the values for A scale with temperature, even though not 
in perfect agreement with the theory. 

We should mention that the s-wave theory we used to fit our data neglects non-linear 
effects that might arise in the vicinity of H c i due to vortex-vortex interaction or due to a 
possible change in the size of the vortex cores. Some experiments performed on s-wave su- 
perconductors have shown non-linear, even \/H, magnetic field dependence for C e \. However 
the field dependence is not consistent for different temperatures, which means depending on 
what temperature the field dependence is probed, the electronic specific heat has a different 
field dependence. Different groups have observed C e i ex H n for almost any value of n be- 
tween 0.5 to 1 depending on what material they studied and at what temperature range they 



performed their experiment f24],|25j . While we can not definitively rule out s-wave symmetry 
as an explanation for our data we believe that a dirty d-wave symmetry gives the most 
consistent and plausible fit to our data. 

In conclusion, our specific heat data strongly suggests that the d-wave symmetry in 
electron doped PCCO at optimal doping is a bulk property of the material. However, due 
to non-magnetic impurities in our sample, the electronic specific heat follows a magnetic 
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field dependence of type Hlog(H) below H*=0.6 T, consistent with dirty d-wave symmetry. 
In addition, the normal state Sommerfeld constant of PCCO, 7„=6.7±0.4 mJ/mole-K 2 , has 
been measured for the first time. 
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FIGURES 

FIG. 1. C/T vs T 2 , where C is the sample total specific heat, at 4 different fields (0, 1, 2, and 
10 Tesla)in the temperature range 2-7 K. The magnetic field is perpendicular to the ab-plane. 

FIG. 2. C(H)-C(0) vs magnetic field, or the field dependent electronic specific heat vs magnetic 
field. A nonlinear behavior is observed below 2T, whereas the high field part has a linear field 
dependence. The electronic specific heat has almost saturated to the normal state value at 8 T. 

FIG. 3. The field dependent electronic specific heat vs magnetic field data at 3.4 K, and the 
theoretical fits to the data. The solid curve is the dirty d-wave fit, the dashed straight line is the 
s-wave fit, and the dashed curve is the clean d-wave fit. 

FIG. 4. The field dependent electronic specific heat vs magnetic field data for 2.4K, 2.7K, and 
3.4K. The lines are the dirty d-wave theory fits to the data. 
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